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SHORTENING IMPULSE RESPONSE FILTER (SIRF) AND DESIGN TECHNIQUE 

THEREFOR 

5 Field of the Invention 

The present invention relates to techniques for shorting the impulse response of 
communication systems, such as discrete multi-tone (DMT) and orthogonal frequency division 
multiplexing (OFDM) communication systems, and more particularly, to methods and apparatus 
for designing a shortening impulse response filter (SIRF). 

10 

Background of the Invention 

It is well known that most communication channels are dispersive in nature and 
introduce a number of distortions. Thus, signals arriving at a receiver are typically corrupted by 
^ intersymbol interference (ISI), crosstalk, echo, and other noise. Thus, receivers must jointly 
45 equalize the channel, to compensate for such intersymbol interference and other distortions, and 
decode the encoded signals at increasingly high clock rates. 

To overcome the effects of intersymbol interference, any two adjacent symbols in 
a conventional DMT or OFDM communication system are separated by a guard period (i.e., a 
3 cyclic prefix (CP)). In addition to providing a mechanism for frame synchronization, the guard 
:?0 interval insures that samples from one symbol block do not interfere with the samples of another 

* block. The length of the impulse response of the physical channel determines the required length 

* of the guard interval. Using a long guard interval, however, reduces the effective throughput of 
the transceiver. Thus, to avoid using a long guard interval, filters are employed to shorten the 
channel impulse response and thereby allowing the use of a shorter guard interval. More 

25 specifically, time domain linear filters, often referred to as shortening impulse response filters 
(SIRFs) or time domain equalizers (TDQs), are typically employed to shorten the channel 
impulse response. 

A number of techniques have been proposed or suggested for designing TDQ 
filters. For a detailed discussion of a number of such filter design techniques, see, for example, 
30 J.W.P. Melsa and R.C. Younce, "Impulse Response Shortening for Discrete Multitone 
Tranceivers," IEEE Trans., COM-44, (12), 1662-1672 (1996); or N. Al-Dahir and J.M. Cioffi, 
"Stable Pole-Zero Modeling of Long FIR Filters With Application to the MMSE-DFE," IEEE 
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Trans., COM-45, (5) 508-513 (1997), each incorporated by reference herein. Generally, these 
filter design algorithms are typically based on least mean square (LMS) or eigenvector calculus. 
While these filter design algorithms are capable of producing very good TDQ filters, they suffer 
from a number of limitations, which if overcome, could greatly improve their ability to shorten 
the channel impulse response and otherwise improve system performance. Specifically, since 
these filter design algorithms have little, if any, control over the frequency response, they may 
produce a frequency response with nulls in the pass-band that degrade the signal-to-noise ratio 
(SNR) of the received signal, translating into a lower bit rate throughput. It has been found, 
however, that removing the nulls in the pass-band is a difficult problem, often requiring a trial 
and error solution. In addition to the null problem, the frequency response in unpredictable and 
severe attenuation and amplification variations could result from call to call. 

A need therefore exists for improved techniques for designing SIRF filters. A 
further need exists for SIRF filters that satisfy constraints in both the time and frequency 
domains to provide improved performance. Yet another need exists for methods and apparatus 
for determining coefficient values for SIRF filters that shorten the channel impulse response in 
the time domain while also providing a frequency response that does not attenuate the received 
signal. 

Summary of the Invention 

Generally, a method and apparatus are disclosed for determining parameters for 
SIRF filters. According to one aspect of the invention, filter coefficients for SIRF filters are 
determined that satisfy constraints in both the time and frequency domains to provide improved 
performance. More specifically, SIRF filters are disclosed that shorten the channel impulse 
response in the time domain while also providing a frequency response that does not attenuate or 
amplify the received signal. 

One or more sets are established to define constraints that the SERF filter must 
satisfy in the time domain, and one or more sets are established to define constraints that the 
SIRF filter must satisfy in the frequency domain. An SIRF filter satisfying both frequency and 
time constraints is obtained by determining the intersection of the various sets. By varying the 
sets utilized to define the time and frequency domain constraints, SIRF filters having a linear or 
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non-linear phase response may be obtained. Vector space projection method is an iterative 
algorithm applied to the sets until the algorithm converges to a solution, i.e., the intersection if 
the sets intersect, or to a point with minimum summing distance to all sets if the sets do not 
intersect. In the case of SIRF design, the solution of the algorithm is the filter coefficients. 
5 A more complete understanding of the present invention, as well as further 

features and advantages of the present invention, will be obtained by reference to the following 
detailed description and drawings. 

Brief Description of the Drawings 

1° FIG. 1 illustrates a variable shortening impulse response filter applied to a signal 

on a dispersive communication channel to shorten the channel impulse response in accordance 
. ^ with the present invention; 

] J FIG. 2 is a schematic block diagram illustrating a modem in which the present 

n invention may be employed. 

t§- 5 FIG. 3 is a flow chart describing the filter coefficient determination process 

] 3 incorporating features of the present invention; 

FIGS. 4A through 4E, collectively, illustrate exemplary pseudo-code for 
V* generating a non-linear SERF filter; and 

% FIG- 5 illustrates the trajectory of iteration in VSPM for two sets until the 

ISO algorithm converges to the intersection, i.e., the solution. 

Detailed Description 

FIG. 1 illustrates the use of a variable shortening impulse response filter 120 on a 
dispersive communication channel 110 to shorten the channel impulse response 130a, b, in 
25 accordance with the present invention. According to one aspect of the present invention, the 
variable SIRF filter 120 can satisfy constraints in both the time and frequency domains to 
provide improved performance. Although described in connection with exemplary DMT and 
OFDM communication systems, it will be understood that the present invention is equally 
applicable to any environment where it is desirable to shorten an impulse response. 
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Vector Space Projection Methods (VSPM) 
According to one aspect of the present invention, vector space projection methods 
are employed to design SIRF filters. For a detailed discussion of vector space projection 
methods (VSPM), see, for example, L. M. Bregman, "Finding the Common Point of Convex Sets 
5 by the Method of Successive Projections, " Dokl. Akad. Nauk USSR, Vol. 162, No. 3, 487 
(1965), incorporated by reference herein. VSPM techniques find a mathematical object (in this 
case a set of coefficients) in a proper vector space that satisfies multiple constraints. When all 
the constraint sets are convex and have a nonempty intersection, the VSPM becomes a powerful 
theory in finding the objects that satisfy all the constraints. As discussed further below, vector 
10 space projection methods employ an iterative algorithm that will converge to a set of finite 
impulse response (FIR) coefficients that satisfies constraints in the time domain (e.g., shortening 
^ the channel impulse response) and constraints in the frequency domain spectrum. 
% Traditionally, VSPM techniques have been employed to design constrained FIR 

3 filters that are tailored to specific applications. See, K.C. Haddad, "Constrained FIR Filter 
|5 Design by the Method of Vector Space Projections," IEEE Trans, on Circuit and Systems II: 
«J Analog and Digital Signal Processing, Vol. 47, No. 8 (Aug. 2000), incorporated by reference 
herein. In the context of the present invention, where VSPM techniques are employed to design 
y an SIRF filter, two (or more) convex sets representing the constraints in time and frequency 
* domains and corresponding projection operators have been mathematically formulated. A first 
2E0 convex set defines the constraints that the SIRF filter 120 must satisfy in the time domain, such 
that when the filter is convolved with the impulse response, the impulse response is shortened. 
Likewise, a second convex set defines the constraints that the SIRF filter 120 must satisfy in the 
frequency domain, such as a low, high or band pass band. p i is defined to be the projection 
operator onto the set q. Thus, to obtain an SIRF filter satisfying both frequency and time 
25 constraints, an intersection of both sets is required. 

Designing the SIRF Filter 
Generally, the present invention designs a variable SIRF filter 120 having an 
impulse response, h, of length # to shorten the impulse response of a channel s , where 

h = (/ 2 (0),/ 2 (l),-,A(iV-l)). 
30 In the frequency domain, h becomes: 
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where A { f0 ) = £ 2 h{n) cos 



n \a> 



and \ a (for linear phase). 



The transformations from the frequency domain to the time domain and vice versa are done 
using the fast Fourier transform and inverse fast Fourier transform, respectively, by discretizing 
co- 

The pass-band and stop-band in the frequency domain are defined to be 
Q p = {co : co p < co < n and Q 5 = {a : 0 < co < co s } respectively. As discussed more fully below, 
the sets involved in designing the SIRF filter 120 in an exemplary embodiment are (i) a time 
domain convex set, Cj, representing the filters with linear phase; (ii) a frequency domain non- 
convex set, C 2 , representing the non-linear phase filters with the appropriate constraints in the 
pass-band and stop-band; (iii) a frequency domain convex set, C 3 , representing the linear phase 
filters with the appropriate constraints in the pass-band and stop-band; (iv) a time domain convex 
set, C 4 , representing all the filters of length jy, and (v) a time domain convex set, C 5 (n), for a 
specific range of values of n , (application dependent). Although C 5 (n) consists of numerous 
convex sets, it is referred to hereinafter as C $ . More specifically, C 5 represents additional 
constraints on the filter h in the time domain. Thus, the time domain sets, C,, C 4 and C 5 , are 
convex, while the frequency domain sets, C 2 and C 3> are convex or non-convex for filters with 
linear or non-linear phase, respectively. The optional frequency domain set, C 3 , constrains the 
filter such that it will have linear phase. 

The sets may be defined mathematically as follows: 

C, = {h e R N :h(n) = h(N-l-n\forn = OX~;N-l\ 

\h eR N :l- a <\H(o>]<l + a for co eQ^ 

C 2 = < 

[and\H(o))<J3 for co eQ s 

hsR N :l-a<A(o))<l + a 
C 3 = \ and <t>(a>) = - co(N - 1)/2 for coeQ p . 
\H(o))<j3 for co eQ s . 
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C 5 (n) = {h eR N :a n < (s *h)„ < p n } (0<n<N+M-l) 
where the vector s referenced in the definition for the set, C 5 , is the impulse response of the 
channel, * denotes convolution, ( s *h)„ denotes the response at time n , and G and p represent 
the desired lower and upper bounds, respectively. M is the size of the discrete channel impulse 
response, s. R N is the Hilbert space of dimension N. 

p i is defined to be the projection operator onto the set q. For a more detailed 
discussion of the computation of the projection operators and the VSPM algorithm generally, 
see, Henry Stark, "Vector Space Projection: A Numerical Approach to Signal and Image 
Processing, Neural Nets, and Optics," (Wiley, 1998), incorporated by reference herein. The two 
iterative algorithms proposed are: 

n * + i = iyV^t, where h 0 is arbitray and projection operators p 2 p 4 and p 5 
are projected onto sets q q q to produce a non-linear phase filter; or 

h k+l =P 1 P 3 P 5 h k , where h 0 is arbitray and projection operators pj p 3 and p 5 
are projected onto sets q q q 5 to produce a linear-phase filter. 

FIG. 2 is a schematic block diagram illustrating a modem 200 in which the 
present invention may be employed. As shown in FIG. 2, the modem 200 includes one or more 
communication ports 230 for receiving a signal from a communication channel 110. As 
previously indicated, the received signal is applied to a variable SIRF filter 120 in accordance 
with the present invention, before being applied to a decoder 240 that decodes the signal in a 
known manner. The coefficients for the SIRF filter 120 are determined by a filter coefficient 
determination process 300, discussed below in conjunction with FIG. 3. 

The filter coefficient determination process 300 may be stored in a data storage 
device 220 that could be implemented as an electrical, magnetic or optical memory, or any 
combination of these or other types of storage devices. Moreover, the term "memory" should be 
construed broadly enough to encompass any information able to be read from or written to an 
address in the addressable space accessed by a processor (not shown). Alternatively, the filter 
coefficient determination process 300 may be embodied on an application specific integrated 
circuit (ASIC). 



Haddad 1 



10 



FIG. 3 is a flow chart describing the filter coefficient determination process 300 
incorporating features of the present invention. As shown in FIG. 3, the filter coefficient 
determination process 300 receives the impulse response of the channel 110 as an input and then 
initializes the SIRF filter 120 to an arbitrary value during step 310 (to provide a starting point). 
Thereafter, the sets that specify the desired filter characteristics are defined during step 320. As 
discussed above, sets q 2 C4 C5 specify the desired characteristics in the time and frequency 
domains for a non-linear filter, while sets q Cj Cj specify the desired characteristics in the 
time and frequency domains for a linear filter. 

The corresponding projection operators p 2 p 4 p 5 or pj p^ p 5 are defined 
during step 350, and are then used to project onto the corresponding sets q C4 Cj or Cj 
C3, Cj during step 360. As shown during step 370, the projections are continued iteratively 
3 until an intersection is reached. The intersection defines the filter coefficients for the SIRF filter 
5; 120, during step 380. Program control then terminates. 

^ FIGS. 4A through 4E, collectively, illustrate exemplary pseudo-code 400 for 

g| 5 generating a non-linear SIRF filter 120. As shown in FIG. 4A, the pseudo-code 400 has an 
J initialization section 410 that initializes a number of parameters and loading the impulse 
^ response for the channel 110. Thereafter, a channel impulse response matrix is established in 
li section 430. The channel impulse response matrix is used for convolution needed for projection 
U onto the set C r As shown in FIG. 4B, the pseudo-code 400 then determines the maximum 
f0 energy of the channel impulse response in section 440. The SIRF filter 120 is initialized to an 

arbitrary value in section 445, and a number of additional parameters are initialized during step 

450. 

As shown in FIG. AC, the first iterative procedure is performed during section 460 
to project onto the set C 2 using the projection operator P 2> The frequency-to-time transformation 

25 is then performed at the end of section 460 using an inverse Fourier transform. As shown in FIG. 
4D, an iterative procedure is performed during section 470 to project onto the set C 4 using the 
projection operator P 4 . As shown in FIGS. 4D and 4E, an iterative procedure is performed 
during section 475 (comprised of sections 475-1 (FIG. 4D) and 475-2 (FIG. 4E) to project the 
projection operator P 5 onto the set C 5 defining additional time characteristics. The time-to- 

30 frequency transformation is then performed at the end of section 475-2 using a Fourier transform. 
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Finally, the determined filter coefficients are applied to the SIRF filter 120 during section 480 
(FIG. 4E). 

FIG. 5 illustrates the trajectory of iteration in VSPM for two exemplary sets, C 1 
and C 2 , until the two sets converge to an intersecting set satisfying the constraints of both sets. 
The solution set C s is the intersection region and x 0 is an arbitrary starting point from which the 
first set is projected onto the second set (at a point defined by P x„). Thereafter, the second set is 
projected onto the first set at a point x 13 where x, equals P 2 P,x 0 . 

It is to be understood that the embodiments and variations shown and described 
herein are merely illustrative of the principles of this invention and that various modifications 
may he implemented by those skilled in the art without departing from the scope and spirit of the 
invention. 



